Hackflex library preparation enables low-cost metagenomic profiling.

Shotgun metagenomic sequencing provides valuable insights into microbial communities, but the high cost of library preparation with standard kits and protocols is a barrier for many. New methods such as Hackflex use diluted commercially available reagents to greatly reduce library preparation costs. However, these methods have not been systematically validated for metagenomic sequencing. Here, we evaluate Hackflex performance by sequencing metagenomic libraries from known mock communities as well as mouse fecal samples prepared by Hackflex, Illumina DNA Prep, and Illumina TruSeq methods. Hackflex successfully recovered all members of the Zymo mock community, performing best for samples with DNA concentrations <1 ng/uL. Furthermore, Hackflex was able to delineate microbiota of individual inbred mice from the same breeding stock at the same mouse facility, and statistical modeling indicated that mouse ID explained a greater fraction of the variance in metagenomic composition than did library preparation method. These results show that Hackflex is suitable for generating inventories of bacterial communities through metagenomic sequencing.

Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures.

BACKGROUND: Social behavior and social organization have major influences on individual health and fitness. Yet, biomedical research focuses on studying a few genotypes under impoverished social conditions. Understanding how lab conditions have modified social organizations of model organisms, such as lab mice, relative to natural populations is a missing link between socioecology and biomedical science. RESULTS: Using a common garden design, we describe the formation of social structure in the well-studied laboratory mouse strain, C57BL/6J, in replicated mixed-sex populations over 10-day trials compared to control trials with wild-derived outbred house mice in outdoor field enclosures. We focus on three key features of mouse social systems: (i) territory establishment in males, (ii) female social relationships, and (iii) the social networks formed by the populations. Male territorial behaviors were similar but muted in C57 compared to wild-derived mice. Female C57 sharply differed from wild-derived females, showing little social bias toward cage mates and exploring substantially more of the enclosures compared to all other groups. Female behavior consistently generated denser social networks in C57 than in wild-derived mice. CONCLUSIONS: C57 and wild-derived mice individually vary in their social and spatial behaviors which scale to shape overall social organization. The repeatable societies formed under field conditions highlights opportunities to experimentally study the interplay between society and individual biology using model organisms.

Assessing co-diversification in host-associated microbiomes.

When lineages of hosts and microbial symbionts engage in intimate interactions over evolutionary timescales, they can diversify in parallel (i.e., co-diversify), producing associations between the lineages' phylogenetic histories. Tests for co-diversification of individual microbial lineages and their hosts have been developed previously, and these have been applied to discover ancient symbioses in diverse branches of the tree of life. However, most host-microbe relationships are not binary but multipartite, in that a single host-associated microbiota can contain many microbial lineages, generating challenges for assessing co-diversification. Here, we review recent evidence for co-diversification in complex microbiota, highlight the limitations of prior studies, and outline a hypothesis testing approach designed to overcome some of these limitations. We advocate for the use of microbiota-wide scans for co-diversifying symbiont lineages and discuss tools developed for this purpose. Tests for co-diversification for simple host symbiont systems can be extended to entire phylogenies of microbial lineages (e.g., metagenome-assembled or isolate genomes, amplicon sequence variants) sampled from host clades, thereby providing a means for identifying co-diversifying symbionts present within complex microbiota. The relative ages of symbiont clades can corroborate co-diversification, and multi-level permutation tests can account for multiple comparisons and phylogenetic non-independence introduced by repeated sampling of host species. Discovering co-diversifying lineages will generate powerful opportunities for interrogating the molecular evolution and lineage turnover of ancestral, host-species specific symbionts within host-associated microbiota.

When microbes go missing: Understanding the impact of diversity loss within the gut microbiome.

Two recent papers published in Cell highlight the power of both top-down and bottom-up approaches to understanding the gut microbiome. The first uses ultra-deep sequencing to identify patterns across a gradient of human industrialization, while the second uses synthetic communities to determine how strain interactions impact microbiome structure and function.

Widespread extinctions of co-diversified primate gut bacterial symbionts from humans.

Humans and other primates harbour complex gut bacterial communities that influence health and disease, but the evolutionary histories of these symbioses remain unclear. This is partly due to limited information about the microbiota of ancestral primates. Here, using phylogenetic analyses of metagenome-assembled genomes (MAGs), we show that hundreds of gut bacterial clades diversified in parallel (that is, co-diversified) with primate species over millions of years, but that humans have experienced widespread losses of these ancestral symbionts. Analyses of 9,460 human and non-human primate MAGs, including newly generated MAGs from chimpanzees and bonobos, revealed significant co-diversification within ten gut bacterial phyla, including Firmicutes, Actinobacteriota and Bacteroidota. Strikingly, ~44% of the co-diversifying clades detected in African apes were absent from available metagenomic data from humans and ~54% were absent from industrialized human populations. In contrast, only ~3% of non-co-diversifying clades detected in African apes were absent from humans. Co-diversifying clades present in both humans and chimpanzees displayed consistent genomic signatures of natural selection between the two host species but differed in functional content from co-diversifying clades lost from humans, consistent with selection against certain functions. This study discovers host-species-specific bacterial symbionts that predate hominid diversification, many of which have undergone accelerated extinctions from human populations.

Home-site advantage for host species-specific gut microbiota.

Mammalian species harbor compositionally distinct gut microbial communities, but the mechanisms that maintain specificity of symbionts to host species remain unclear. Here, we show that natural selection within house mice (Mus musculus domesticus) drives deterministic assembly of the house-mouse gut microbiota from mixtures of native and non-native microbiotas. Competing microbiotas from wild-derived lines of house mice and other mouse species (Mus and Peromyscus spp.) within germ-free wild-type (WT) and Rag1-knockout (Rag1-/-) house mice revealed widespread fitness advantages for native gut bacteria. Native bacterial lineages significantly outcompeted non-native lineages in both WT and Rag1-/- mice, indicating home-site advantage for native microbiota independent of host adaptive immunity. However, a minority of native Bacteriodetes and Firmicutes favored by selection in WT hosts were not favored or disfavored in Rag1-/- hosts, indicating that Rag1 mediates fitness advantages of these strains. This study demonstrates home-site advantage for native gut bacteria, consistent with local adaptation of gut microbiota to their mammalian species.

Convergent remodelling of the gut microbiome is associated with host energetic condition over long-distance migration.

The gut microbiome can be thought of as a virtual organ given its immense metabolic capacity and profound effects on host physiology. Migratory birds are capable of adaptively modulating many aspects of their physiology to facilitate long-distance movements, raising the hypothesis that their microbiome may undergo a parallel remodeling process that helps to meet the energetic demands of migration.To test this hypothesis, we investigated changes in gut microbiome composition and function over the fall migration of the Blackpoll Warbler (Setophaga striata), which exhibits one of the longest known autumnal migratory routes of any songbird and rapidly undergoes extensive physiological remodeling during migration.Overall, our results showed that the Blackpoll Warbler microbiome differed significantly across phases of fall migration. This pattern was driven by a dramatic increase in the relative abundance of Proteobacteria, and more specifically a single 16S rRNA gene amplicon sequence variant belonging to the family Enterobacteriaceae. Further, Blackpoll Warblers exhibited a progressive reduction in microbiome diversity and within-group variance over migration, indicating convergence of microbiome composition among individuals during long-distance migration. Metagenomic analysis revealed that the gut microbiome of staging individuals was enriched in bacterial pathways involved in vitamin, amino acid, and fatty acid biosynthesis, as well as carbohydrate metabolism, and that these pathways were in turn positively associated with host body mass and subcutaneous fat deposits.Together, these results provide evidence that the gut microbiome of migratory birds may undergo adaptive remodeling to meet the physiological and energetic demands of long-distance migration.

Influential factors of saliva microbiota composition.

The oral microbiota is emerging as an influential factor of host physiology and disease state. Factors influencing oral microbiota composition have not been well characterised. In particular, there is a lack of population-based studies. We undertook a large hypothesis-free study of the saliva microbiota, considering potential influential factors of host health (frailty; diet; periodontal disease), demographics (age; sex; BMI) and sample processing (storage time), in a sample (n = 679) of the TwinsUK cohort of adult twins. Alpha and beta diversity of the saliva microbiota was associated most strongly with frailty (alpha diversity: ß = -0.16, Q = 0.003, Observed; ß = -0.16, Q = 0.002, Shannon; ß = -0.16, Q = 0.003, Simpson; Beta diversity: Q = 0.002, Bray Curtis dissimilarity) and age (alpha diversity: ß = 0.15, Q = 0.006, Shannon; ß = 0.12, Q = 0.003, Simpson; beta diversity: Q = 0.002, Bray Curtis dissimilarity; Q = 0.032, Weighted UniFrac) in multivariate models including age, frailty, sex, BMI, frailty and diet, and adjustment for multiple testing. Those with a more advanced age were more likely to be dissimilar in the saliva microbiota composition than younger participants (P = 5.125e-06, ANOVA). In subsample analyses, including consideration of periodontal disease (total n = 138, periodontal disease n = 66), the association with frailty remained for alpha diversity (Q = 0.002, Observed ASVs; Q = 0.04 Shannon Index), but not beta diversity, whilst age was not demonstrated to associate with alpha or beta diversity in this subsample, potentially due to insufficient statistical power. Length of time that samples were stored prior to sequencing was associated with beta diversity (Q = 0.002, Bray Curtis dissimilarity). Six bacterial taxa were associated with age after adjustment for frailty and diet. Of the factors studied, frailty and age emerged as the most influential with regards to saliva microbiota composition. Whilst age and frailty are correlates, the associations were independent of each other, giving precedence to both biological and chronological ageing as processes of potential importance when considering saliva microbiota composition.

Microbiomes: Infant Chimps Crawling with Bacteria.

Human guts are colonized at birth by a limited set of microbes that gradually increases in diversity throughout infancy. A new study reports the opposite pattern in infant chimpanzees, raising questions about how host-microbiota relationships have changed during ape evolution.

Microbiota assembly, structure, and dynamics among Tsimane horticulturalists of the Bolivian Amazon.

Selective and neutral forces shape human microbiota assembly in early life. The Tsimane are an indigenous Bolivian population with infant care-associated behaviors predicted to increase mother-infant microbial dispersal. Here, we characterize microbial community assembly in 47 infant-mother pairs from six Tsimane villages, using 16S rRNA gene amplicon sequencing of longitudinal stool and tongue swab samples. We find that infant consumption of dairy products, vegetables, and chicha (a fermented drink inoculated with oral microbes) is associated with stool microbiota composition. In stool and tongue samples, microbes shared between mothers and infants are more abundant than non-shared microbes. Using a neutral model of community assembly, we find that neutral processes alone explain the prevalence of 79% of infant-colonizing microbes, but explain microbial prevalence less well in adults from river villages with more regular access to markets. Our results underscore the importance of neutral forces during microbiota assembly. Changing lifestyle factors may alter traditional modes of microbiota assembly by decreasing the role of neutral processes.

Treatment-Specific Composition of the Gut Microbiota Is Associated With Disease Remission in a Pediatric Crohn's Disease Cohort.

BACKGROUND: The beneficial effects of antibiotics in Crohn's disease (CD) depend in part on the gut microbiota but are inadequately understood. We investigated the impact of metronidazole (MET) and metronidazole plus azithromycin (MET+AZ) on the microbiota in pediatric CD and the use of microbiota features as classifiers or predictors of disease remission. METHODS: 16S rRNA-based microbiota profiling was performed on stool samples from 67 patients in a multinational, randomized, controlled, longitudinal, 12-week trial of MET vs MET+AZ in children with mild to moderate CD. Profiles were analyzed together with disease activity, and then used to construct random forest models to classify remission or predict treatment response. RESULTS: Both MET and MET+AZ significantly decreased diversity of the microbiota and caused large treatment-specific shifts in microbiota structure at week 4. Disease remission was associated with a treatment-specific microbiota configuration. Random forest models constructed from microbiota profiles before and during antibiotic treatment with metronidazole accurately classified disease remission in this treatment group (area under the curve [AUC], 0.879; 95% confidence interval, 0.683-0.9877; sensitivity, 0.7778; specificity, 1.000; P < 0.001). A random forest model trained on pre-antibiotic microbiota profiles predicted disease remission at week 4 with modest accuracy (AUC, 0.8; P = 0.24). CONCLUSIONS: MET and MET+AZ antibiotic regimens in pediatric CD lead to distinct gut microbiota structures at remission. It may be possible to classify and predict remission based in part on microbiota profiles, but larger cohorts will be needed to realize this goal.

Role of priority effects in the early-life assembly of the gut microbiota.

Understanding how microbial communities develop is essential for predicting and directing their future states. Ecological theory suggests that community development is often influenced by priority effects, in which the order and timing of species arrival determine how species affect one another. Priority effects can have long-lasting consequences, particularly if species arrival history varies during the early stage of community development, but their importance to the human gut microbiota and host health remains largely unknown. Here, we explore how priority effects might influence microbial communities in the gastrointestinal tract during early childhood and how the strength of priority effects can be estimated from the composition of the microbial species pool. We also discuss factors that alter microbial transmission, such as delivery mode, diet and parenting behaviours such as breastfeeding, which can influence the likelihood of priority effects. An improved knowledge of priority effects has the potential to inform microorganism-based therapies, such as prebiotics and probiotics, which are aimed at guiding the microbiota towards a healthy state.

Characterization of the facial microbiome in twins discordant for rosacea.

Previously, we determined that genetic and environmental factors contributed equally towards rosacea in twins. To assess an environmental factor, we characterized the malar cheek bacterial microbiome from twins discordant for rosacea. We found no significant difference in facial microbiome alpha and beta diversity between related twins discordant for rosacea. However, the relative percentage abundance of Gordonia and Geobacillus, low-abundant genera, was positively and negatively associated with rosacea severity, respectively. Our data demonstrate a significant correlation between facial microbiome and severity of rosacea in genetically matched twins and importantly that overall microbiome composition is largely unchanged.

Comparison of pectin-degrading fungal communities in temperate forests using glycosyl hydrolase family 28 pectinase primers targeting Ascomycete fungi.

Fungi have developed a wide assortment of enzymes to break down pectin, a prevalent polymer in plant cell walls that is important in plant defense and structure. One enzyme family used to degrade pectin is the glycosyl hydrolase family 28 (GH28). In this study we developed primers for the amplification of GH28 coding genes from a database of 293 GH28 sequences from 40 fungal genomes. The primers were used to successfully amplify GH28 pectinases from all Ascomycota cultures tested, but only three out of seven Basidiomycota cultures. In addition, we further tested the primers in PCRs on metagenomic DNA extracted from senesced tree leaves from different forest ecosystems, followed by cloning and sequencing. Taxonomic specificity for Ascomycota GH28 genes was tested by comparing GH28 composition in leaves to internal transcribed spacer (ITS) amplicon composition using pyrosequencing. All sequences obtained from GH28 primers were classified as Ascomycota; in contrast, ITS sequences indicated that fungal communities were up to 39% Basidiomycetes. Analysis of leaf samples indicated that both forest stand and ecosystem type were important in structuring fungal communities. However, site played the prominent role in explaining GH28 composition, whereas ecosystem type was more important for ITS composition, indicating possible genetic drift between populations of fungi. Overall, these primers will have utility in understanding relationships between fungal community composition and ecosystem processes, as well as detection of potentially pathogenic Ascomycetes.

The effect of microbial colonization on the host proteome varies by gastrointestinal location.

Endogenous intestinal microbiota have wide-ranging and largely uncharacterized effects on host physiology. Here, we used reverse-phase liquid chromatography-coupled tandem mass spectrometry to define the mouse intestinal proteome in the stomach, jejunum, ileum, cecum and proximal colon under three colonization states: germ-free (GF), monocolonized with Bacteroides thetaiotaomicron and conventionally raised (CR). Our analysis revealed distinct proteomic abundance profiles along the gastrointestinal (GI) tract. Unsupervised clustering showed that host protein abundance primarily depended on GI location rather than colonization state and specific proteins and functions that defined these locations were identified by random forest classifications. K-means clustering of protein abundance across locations revealed substantial differences in host protein production between CR mice relative to GF and monocolonized mice. Finally, comparison with fecal proteomic data sets suggested that the identities of stool proteins are not biased to any region of the GI tract, but are substantially impacted by the microbiota in the distal colon.

Use of 16S rRNA sequencing and quantitative PCR to correlate venous leg ulcer bacterial bioburden dynamics with wound expansion, antibiotic therapy, and healing.

Clinical diagnosis of infection in chronic wounds is currently limited to subjective clinical signs and culture-based methods that underestimate the complexity of wound microbial bioburden as revealed by DNA-based microbial identification methods. Here, we use 16S rRNA next generation sequencing and quantitative polymerase chain reaction to characterize weekly changes in bacterial load, community structure, and diversity associated with a chronic venous leg ulcer over the 15-week course of treatment and healing. Our DNA-based methods and detailed sampling scheme reveal that the bacterial bioburden of the wound is unexpectedly dynamic, including changes in the bacterial load and community structure that correlate with wound expansion, antibiotic therapy, and healing. We demonstrate that these multidimensional changes in bacterial bioburden can be summarized using swabs taken prior to debridement, and therefore, can be more easily collected serially than debridement or biopsy samples. Overall, this case illustrates the importance of detailed clinical indicators and longitudinal sampling to determine the pathogenic significance of chronic wound microbial dynamics and guide best use of antimicrobials for improvement of healing outcomes.

Evolutionary analysis of glycosyl hydrolase family 28 (GH28) suggests lineage-specific expansions in necrotrophic fungal pathogens.

Glycosyl hydrolase family 28 (GH28) is a set of structurally related enzymes that hydrolyze glycosidic bonds in pectin, and are important extracellular enzymes for both pathogenic and saprotrophic fungi. Yet, very little is understood about the evolutionary forces driving the diversification of GH28s in fungal genomes. We reconstructed the evolutionary history of family GH28 in fungi by examining the distribution of GH28 copy number across the phylogeny of fungi, and by reconstructing the phylogeny of GH28 genes. We also examined the relationship between lineage-specific GH28 expansions and fungal ecological strategy, testing the hypothesis that GH28 evolution in fungi is driven by ecological strategy (pathogenic vs. non-pathogenic) and pathogenic niche (necrotrophic vs. biotrophic). Our results showed that GH28 phylogeny of Ascomycota and Basidiomycota sequences was structured by specific biochemical function, with endo-polygalacturonases and endo-rhamnogalacturonases forming distinct, apparently ancient clades, while exo-polygalacturonases are more widely distributed. In contrast, Mucoromycotina and Stramenopile sequences formed taxonomically-distinct clades. Large, lineage-specific variation in GH28 copy number indicates that the evolution of this gene family is consistent with the birth-and-death model of gene family evolution, where diversity of GH28 loci within genomes was generated through multiple rounds of gene duplication followed by functional diversification and loss of some gene family members. Although GH28 copy number was correlated with genome size, our findings suggest that ecological strategy also plays an important role in determining the GH28 repertoire of fungi. Both necrotrophic and biotrophic fungi have larger genomes than non-pathogens, yet only necrotrophs possess more GH28 enzymes than non-pathogens. Hence, lineage-specific GH28 expansion is the result of both variation in genome size across fungal species and diversifying selection within the necrotrophic plant pathogen ecological niche. GH28 evolution among necrotrophs has likely been driven by a co-evolutionary arms race with plants, whereas the need to avoid plant immune responses has resulted in purifying selection within biotrophic fungi.

Circadian input kinases and their homologs in cyanobacteria: evolutionary constraints versus architectural diversification.

The circadian input kinase A (cikA) gene encodes a protein relaying environmental signal to the central circadian oscillator in cyanobacteria. The CikA protein has a variable architecture and usually consists of four tandemly arrayed domains: GAF, histidine kinase (HisKA), histidine kinase-like ATPase (HATPase_c), and a pseudo-receiver (REC). Among them, HisKA and HATPase_c are the least polymorphic, and REC is not present in heterocystic filamentous cyanobacteria. CikA contains several conserved motifs that are likely important for circadian function. There are at least three types of circadian systems, each of which possesses a different set of circadian genes. The originally described circadian system (kaiABC system) possesses both cikA and kaiA, while the others lack either only cikA (kaiABC (Delta)) or both (kaiBC). The results we obtained allowed us to approximate the time of the cikA origin to be about 2600-2200 MYA and the time of its loss in the species with the kaiABC (Delta) or kaiBC system between 1100 and 600 MYA. Circadian specialization of CikA, as opposed to its non-circadian homologs, is a result of several factors, including the unique conserved domain architecture and high evolutionary constraints of some domains and regions, which were previously identified as critical for the circadian function of the gene.